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We study a new minimal scalar-tensor model of gravity with Brans-Dicke factor 
cj(<J>) = and cosmological factor n($). The constraints on n($) from known 
gravitational experiments are derived. We show that almost any time evolution of 
the scale factor in a homogeneous isotropic Universe can be obtained via a properly 
chosen n($) and discuss the general properties of models of this type. 
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At present general relativity (GR) is the most successful theory of gravity in describinggravi- 
tational phenomena at laboratory-, earth-surface-, solar-system- and star-systems scales Al- 
^ though it gives quite good a description of these phenomena at galaxies scales and at the scales 
of the whole Universe, there exist well known problems such as the rotation of galaxies, the initial 
^ ' singularity problem, the early Universe, the recent discovery of the accelerated expansion of the 
Universe, etc. It is likely, that the indicated problems call for an extension of the GR framework. 
The most promising modern theories which incorporate naturally GR like supergravity and 
, (super)string theories, unfortunately introduce, at least at their present stage of development, a 
large number of new fields without any real physical basis and are far from being experimentally 
testable. Therefore, it seems sensible to explore some other models which are (i) compatible with 
the established gravitational experiments, (ii) obtained through a minimal extension of GR, (iii) 
^T) [ promising candidates for overcoming the above problems, and which (iv) may be considered as a 
part of some more general modern theories. In the present article we outline the general properties 
of one such model which is built upon an additional scalar field and that differs from the known 
, inflationary and quintessential models. 

■ We call minimal dilatonic gravity (MDG) the scalar-tensor model of gravity with an action 

Ag,a = -^I d*xy^mR + 2AU{^)). (1) 

bJT), Here, the well known Brans-Dicke factor has been denoted by uj{^) = 0, A stays for the cosmological 
constant, while n(<I>) is the diniensionless cosmological factor. The matter action. Am, and the 
• , matter equations of motion will have the usual GR form and do not dependend directly on the 
^ ■ dilaton field $. 

^ , The field equations for the metric gap and the dilaton field $ are given by: 

$(G„^-An($)g„;3)~(VaV;3-g„^n)$ = ^r„^, 

□$+A||($)=5f^r (2) 

and yield the energy-momentum conservation law, Vq Tp = 0. 

The variation of the action (Q) with respect to the metric gaff lieds to the first one of the equations 
(^. The variation of ([^) with respect to the dilaton field <I> produces the following local algebraic 
relation between this field and the scalar curvature R: 

i?-h2A($^($) +n($)) = 0. (3) 

In combining the latter relation with the trace of the first of the equations (||), the field equation 
for the dilaton $ is obtained. Because of the absence of the standard kinetic term ~ aj(<I>)(V$)^ 
for the dilaton in action (|^), its dependence on nonzero space-time curvature (R^O) is essential 
for the existence of this equation. There, ^ — |$($^ — 11) is the derivative of the dilatonic 
potential V{^). 
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The presence of a potential term for the scalar fields in the action of the theory, was considered 
many times from different points of view. For example, the extended inflation models and the 
quintessential models contain such a term, see ^ and the references therein. 

In Brans-Dicke's theory, n($) = 0, and one has to make a choice for the function ijj(<i>) only, 
in order to fix the model. However, recent astrophysical data from type la supcrnovac at high 
red-shifts, on the one hand, as well as from cosmic microwave background, on the other hand, lean 
a strong support to the presence of an additional A-term in the action of the theory Q . 

In more general scalar-tensor theories of gravity, two unknown functions, uj{^) and n(<I>), are to 
be fixed. The old attempts to extract from the known gravitational experiments simultaneously 
with the dilaton mass ^ also a constant parameter u = const ^ 0, turned out to be less 
successful, living great uncertainties for the corresponding values of the above parameters, see 
A new promising method to determine a;($) and n($) (or some equivalent functions) using 
astrophysical observations, was developed in the recent work |^]. Unfortunately, at present, we 
don't have sufficiently reliable data, as required for applying this method, so that, at present, the 
form of both functions still remains an open problem. 

As a scalar-tensor theory with only one unknown function n(<l>), the MDG is much better defined. 
Its investigation was pioneered by O'Hanlon in connection with Fujii's theory of massive dilatons 
as early as in |7|, however, without any connection with the cosmological constant problem. 

In general, such a minimal model may be related with the so called nonlinear theories of gravity 
(see [|| and the references therein). Indeed, the variation of the action (|l|) with respect to the dilaton 
field $ gives the local algebraic relation (H) between scalar curvature R and dilaton $ in place of the 
differ ential equation, which is due to the absence of the Brans-Dicke kinetic term ~ aj(<I')(V$)^ 
in the action for = 0. If ($n($)) ^0, one can solve the relation (^ with respect to the 
field $ and it becomes a local function of the scalar curvature: <I> = $(i?). (For w ^ the last 
relation will have non local integral form.) The substitution of this function back into the action 
(|l|) transforms it to the action of nonlinear theories of gravity Ang — S d'^^^/Wi fi^) with 
f{R) = i?$(i?) + 2A$(i?)n($(i?)). Hence, the action (|l|) is of the Hclmhohz type for nonlinear 
gravity g. 

MDG may be considered also as a part of more general theories of gravity like the (4-1-1)- 
dimensional Kaluza-Klein model described by Fujii in or, the model of gravity with violated 
local conformal symmetry in affine-connected spaces, which probably may be related to string 
theory @. 

In this special scalar-tensor model of gravity the cosmological factor H($) is the only function 
which has to be chosen to comply with gravitational experiments and astrophysical observations. 
In the present article, we shall demonstrate that this results in much better defined predictions for 
the dilaton mass, which is now derived from known gravitational experiments. It turns out that 
the results of these experiments on Earth surface and in solar, or star systems depend mainly on 
the mass term in the dilaton potential. A numerical analysis of the boson stars in massive dilatonic 
gravity based on MDG was performed very recently It shows that the boson star structure 
is sensitive, typically within a few percent accuracy only, to the dilaton mass rather than to the 
exact form of the potential. Hence, we need some other source of information about the form of 
the cosmological factor. 

In the present article we shall show how one can solve the inverse cosmological problem in 
MDG. By inverse cosmological problem we mean the determination of the cosmological factor 
H(<I>) which yields a given time evolution of the scale parameter A[t) in the Robertson- Walker 
(RW) model of the Universe. Such an approach for the determination of the cosmological factor 
may be considered as a further development of the general idea to find the dilaton potential from 
astrophysical observations as discussed first in the articles There, the relation between the 

dilaton potential and the luminosity distance was obtained from observations of type la supernovae 
explosions at high red-shifts. Similar consideration based on real observational data may take place 
in the model of MDG, too. Our approach, in being based upon the evolution of the cosmological 
factor, gives another possibility to extract the dilaton potential from data and seems to be more 
profound. 

An essential new element of our MDG is the nonzero constant A which we identify with standard 
cosmological constant. The new astrophysical data J^a = -65 ± .13, Hq = (65 ± 5) km s~^Mps~^ Q 
give observed value A"*"* — 3^aHqC~'^ = (.98 ± .34) x 10^*^^0771^^. Though the present confidence 
level of SIa is not as high, we accept its value as a basic quantity which defines the natural units 
for all cosmological quantities: the cosmological length Ac — l/VA"'"* = (1.02 ± .18) x 10^^c77i, the 
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cosmological time Tc := A^jc = (3.4±.6) x lO^^s = (10.8±1.9)Gt/r, the cosmological energy density 
^ = (1.16±.41) X cto"^ the cosmological energy := SA^Ec = BA-^/^c^k-^ = 

(3.7 ± .7) X the cosmological momentum Pc := 3c/(k\/A°^) = (1.2 ± .2) x lO^^g cm s^^ 

and new the cosmological unit for the action Ac ■— 3c/(kA°^*) = (1.2 ± .4) x 10^^^ h, k being 
Einstein constant. Further, we use dimensionless variables like; t :~ t/Tc, a :— A/Ac, h :— HTc 
(H :— A^^dA/dt being Hubble parameter), Cc = ec/|ec| = ±1, e := e/jed-matter energy density, 
etc. 



II. SOLAR SYSTEM AND EARTH-SURFACE GRAVITATIONAL EXPERIMENTS 

We use the well studied scalar-tensor theories of gravity ||l[] , Q , [|j , to derive the properties 
of cosmological factor n(<i>) and predictions for gravitational experiments. 

A. General Consideration 

1. The MDG with A = is a Brans-Dicke theory with oj — and contradicts the data from 
the solar system gravitational experiments, because the value ^ of the post Newtonian parameter 
7 = is far from the experimental constraint 7 = 1 ± .001 Q|. We must introduce cosmological 
term An($) 7^ in action (|l|) to overcome this problem. 

2. In contrast to O'Hanlon's model we wish MDG to reproduce GR with A ^ for <i> = 
$ = const 0. Then from action (|l|) we obtain normalization condition n($) — 1 and Einstein 
constant n = R/^. 

3. In vacuum, far from matter MDG have to allow weak field approximation <i> — $(1+C), |C| ^ 1 
which we consider in harmonic gauge. Then the second of equations (|^), written in linearized form 
with respect to C: 

gives ^($) = 1/$. Taylor series expansion of the function ^($) around the value $ starts with 
term which can be connected with the dilaton mass m$ = Z$ being the dilaton Compton 
length. Thus we obtain ^(l>) = |p^2^-2 ^^^^ 

n=i+C+^C'+o(c=^), (4) 

p — ^ being a new dimensionless parameter. 

4. For few point particles of masses Too being source of metric and dilaton fields in equations 
(^ we obtain Newtonian approximation ((?a/3 = ^a/3+/ia/3i |^a/3|^l) for the gravitational potential 
93 (r) and the dilaton field <i>(r): 

^(r)/c2 = -i|^^(l + a(p)e-l-'-'-l/'-) 

a 

a 

$(r)/$ = l+2 G^y^s^e-I"'-'"-!/'*, (5) 

^ a 

G^^{l — ^p^) is Newton's constant, M — J2a The term 

a a 

in (p{r) is known from GR with A 7^ 0. It represents an universal anty-gravitational interaction of 
test mass to with any mass nia via repellent elastic force 

F,o = iATOc2^(r-ro). (6) 
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For solar system distances / < 1000 AU the whole second term in ip may be neglected ||T^ , being 
of order < 10~^^. Then we arrive to the known form of gravitational potential (p but with 
specific for MDG coefiicient 

«b) = 

and the comparison of both possibilities: a > |, or a < — 1 with Cavendish type experiments yields 
the experimental constraint < 1.6 [mm] {p < 2 x 10^^^), see the articles by De Riijula and by 
Fischbach and Talmadge in Q and the references therein. Hence, in the solar system phenomena 
the factor e~'/'* has fantastically small values (< exp{—10^^) for the Earth-Sun distances I, or 
< exp{—3 X 10^°) for the Earth-Moon distances /) and there is no hope to find some differences 
between MDG and GR in this domain. 

The constraint m$c^ > 10^'' [ey] does not exclude a small value (with respect to the elementary 
particles scales) for the rest energy of the hypothetical ^-particle. 

5. The parameterized-post-Newtonian (PPN) solution of equation (||) is complicated, but be- 
cause of the constraint p < 10~^* we can neglect the second term in ip, set a = | and use with great 
precision Helbig's PPN formalism which differs essentially from the standard one 0] for zero mass 
dilaton fields. We obtain much more definite predictions then general relations between a and the 
length 1$ given in the articles by De Riijula and by Helbig 1^, because of the condition uj = 0. 

B. The Basic Gravitational Effects in MDG 

A brief treatment of basic gravitational effects gives: 

1. Nordtvedt Effect 

In MDG a body with significant gravitational self-energy = J^t^c ^ k'^'' ^^^^ ^'-'^ move 
along geodesies due to additional universal anty- gravitational force: 

F« = -l^^oV*. (7) 

For usual bodies it is too small even at distances |r — TqI < because of the small factor Eq- 
Hence, in MDG we have no strict strong equivalence principle. Nevertheless, the week equivalence 
principle is not violated. 

The experimental data for Nordtvedt effect caused by the Sun are formulated as a constraint 
77 = ± .0015 jl] on the parameter 77 which in MDG becomes a function of the distance I to the 
source: 

77(0 - -1(1 + V'*)e-'/'*. 

Taking into account the value of the Astronomical Unit (AU) l^^ « 1.5 x 10^^ [to] we obtain from 
the experimental value of rj the constraint < 2 x 10^° [to]. 

2. Time Delay of Electromagnetic Waves 

The standard action for electromagnetic field and the Maxwell equations do not depend directly 
on the field $. Therefore the influence of this field on the electromagnetic phenomena like motion 
of electromagnetic waves in vacuum is possible only indirectly - via its influence on the space-time 
metric. The solar system measurements of the time delay of the electromagnetic pulses give the 
above used value 7 = 1 ± .001 of this post Newtonian parameter Q]. In MDG this yields the 
relation (1 ± .001)g(Z^„) — 1 and gives once more the constraint Z$ < 2 x 10^''[to]. Here 

g{l) := l + i(l + ///$)e-'/'*. 
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3. Perihelion Shift 

Helbig's results (s) applied in MDG give the following formula for the perihelion shift of a planet 
orbiting around the Sun (with mass Mq): 

Here Ip is the semi-major axis of the orbit of planet and 

fc(/p) « 1 + ^ (4 + e-'^/'- ~ 

is obtained neglecting its eccentricity. The observed value of perihelion shift of Mercury gives the 
constraint l,p < 10^ [m] (see the article by De Riijula ||]). 

The above week restrictions on derived in MDG from gravitational experiments show that in 
presence of dilaton field $ are impossible essential deviations from GR in solar system. Formulae 
show that observable deviations from Newton law of gravity can not be expected at distances 
greater then few cm, too. 



III. APPLICATION OF MDG IN COSMOLOGY 
A. General Equations of MDG for RW Universe 

The design of realistic model of the Universe lies beyond the scope of the present article. Here 
we wish only to outline some general features of the MDG, applied to cosmological problems and 
to show that in general it is able to correct the well known shortcomings of GR in this domain. We 
derive the equations of the inverse cosmological problem in MDG and demonstrate some simple 
solutions, which give indications for unexpected new physics. 

Consider RW adiabatic homogeneous isotropic Universe with dsj^^, = c^dt^ — A^dl^, t = TcT, 
A{t) — Aca^r) and dimensionless dlf. = jz!^ +P{d6'^ + sin'^9)dip'^ {k = —1,0,1) in presence of 
matter with energy-density e = £ce(a)/^ and pressure p — ScPeia)/^. 

From equations (||) and (||) we obtain the following basic dynamical equations of MDG for RW 
Universe (see for details the Appendix): 

1 d'^a _i_ 1 f da 



a dT^ 



+ i7(f)^ + ^ = i(*^(*) + n($)): 



Using Hubble parameter h{a) = a~^^(T(a)) (T(a) is the inverse function of a(T)), new variable 
A = In a and denoting by prime the differentiation with respect to A we write down the equations 
(^ as a second order non- autonomous system for the functions $(A) and h^{\): 

i(/.2)' + 2h^ + = i (ci>^($) + H($)) , 

4- {h^ + ke-^^) $ = i(«>H($) + e(e^)). (9) 

For given function H($) the solutions of system (^ determine the functions <i>(a) and h{a) and 
the time dependence of the scale parameter a(r) is given implicitly by the relation 

r(a) / da/{ah{a)) + (10) 

Instead, we can exclude the cosmological factor H($) arriving to linear differential equation of 
second order for the function <i>(A): 

$" + (^-1) $' + 2 kh-\-^^) 3i,e'. (11) 



In terms of the function ?A(a) ~ yj\h{a)\/a $(a) it reads: 

^" + = 5, (12) 
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where we have introduced the new functions 



~ 2 h A \ h I 2 h ^ h^^ ' 



(13) 

For the important coefhcient n one can obtain the representation 

where —q = = 1 + ^ is decelerating parameter and S'[a(T)] = — I is 

Schwartz derivative of the scale factor a{T). 



B. The Inverse Cosmological Problem 



Now we are ready to consider the inverse cosmological problem: to find the corresponding 
cosmological factor n($) and dilatonic potential V{^) which leads to a given evolution of the 
Universe, determined by function a(T). It is remarkable that in MDG an unique solution of this 
problem exist for almost any three times differentiablc function a(T). Indeed, for given a(T) we 
may construct a function h{X) and find the general solution ^(A, Ci, C2) of the linear second order 



differential equation (11). It depends on two arbitrary constants Ci^2 which have to be determined 
from the additional conditions n(l>) = 1, ^(i") = ^(^) = They ensure the 

self-consistence of the MDG model yielding initial conditions for equation (^ij) at point A which is 
real solution of the algebraic equation 

r(A) = -4. (14) 

This is a new form of relation (|) and r(A) = -6 (K^i^)' + 2h'^ + fce-^^) is the dimensionless scalar 
curvature r — R/A. Then: 

$ = -4e-(l + |p2) /{j'oo{l + lp')+4p'h^r') , 

|.7$ = -|pV/(l + |p2). (15) 

Here and later on, the bar sign means that the corresponding quantities are calculated at the 
point A, and joo = Goo /A = 3 [h^ + ke~'^^^ is the dimensionless 00-component of Einstein tensor. 
Hence, the values of all "bar" quantities (including R in action (|l|)) may be determined from the 
time evolution a(r) of the Universe via the solution A = In a of the equation (^. 

In their turn the quantities $ and determine the values of the constants Ci^2 and an unique 
solution <!>( A) of the equation (pT|). Indeed, let $i(A) and<I>2(A) are fundamental system of solutions 
of the homogeneous equation associated with non-homogeneous one (|ll|). Then A(A) := <I>i<I>2 — 
$2$'i = (AVa)e^ h{\) + and 

$(A) = Ci$i(A) + C2$2(A)-l-$e(A) (16) 

is the general solution of equation ( [Til ) which depends on two constants 

Ci = ($2^ - $2$')/^> C2 = ($1$' - $il>)/A 

and on contribution of matter described by the term: 

de<i>i/{ah) - $1 / 



The cosmological factor H and the potential V as functions of the variable A are determined by 

the equations 

H(A) = Joo + 3/i2$7$ - e/$, v{X) ^ l j ^ (*n' - $'n) dX (17) 

which define the functions H(<f>) and V{^) implicitly, as well. 
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C. Simple Examples 



Let us demonstrate the above procedure for solution of the inverse cosmological problem on 
several simple examples: 

1. Evolution law a(T) = (wr)"'^/'^, (where cj is a free parameter) gives 

h{X)^^e-''^, q = v-l, -n2(A) = iii^f±^ + 2A:5^e2(''-i)^ 

and the equation 

-2 30^-^-2(1-.)^^ 

for a. Then: 

i) For > 2 we have real solution a only if k = +1: 

$i(a) = a'^I^{ba''-') , $2(0) = a^K^{ba''-^) , 



■= 2(i^-'i)"'^ ~ being the order of Bessel functions Ifi,Jp,, K^, and b := \Jj^^^^- For GR-hke 

law a^y/r {v — 2) in MDG we obtain positive value (fc — +1) for the three-space curvature, A = 0, 
/X = I and Bessel functions are reduced to elementary functions, 
ii) When < 2 all values k = —1,0, +1 are admissible: 

- for k = +l the solutions $1.2 are the same as above; 

- for A; = we have 

- for k — —l the solutions are: 

In the special case of linear evolution a^r {v = 1) — n^(A) = 3 + 

$i,2(a)-ai±^ 

and the root A = i In |(fc + o;^) is real for all values of cj^ > 0, if fc = 0, +1. For k = —l the root A 
will be real if \u)\ > 1. 

2. Evolution law a(T) = e'^^ gives 

h{X)=cj, g = -n2(A)-i + 



and the equation |a (1 — Scj ) ~ k whith root a= \J 2\i-3i^'^\ ■ Now we have the following solutions: 

i) For |w|<^ and k = +\: 

*i=«cosh(^), cl>2 = asinh(^) ; 

ii) For |t^|>^ and fc = -l: 

cos(|^) , $2 = a sin(|^) • 

Conditions (|lj) and (|l5| ) exclude exact exponential expansion of spatially flat Universe (fc = ). 
For inflationary scenario in this case one may use a scale factor a{T) = e^'^ + const which turns out 
to be possible if const ^Q. 
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IV. SOME SPECIFIC PROPERTIES AND POSSIBLE FURTHER DEVELOPMENTS OF 

THE MDG MODEL 



In the following concluding remarks we stress some basic properties of MDG models and outline 
some possible important further developments to be considered elsewhere. First of all we have 
to stress that actually here we consider the general properties of certain class of such models. 
Depending on the form of cosmological factor n($) the concrete representatives of this class of 
models of extended gravity will have different physical properties and consequences, having in the 
same time some general features. 

As we see, the MDG models of the RW Universe have the following specific properties: 

1) In domains where n>0 the dilatonic field oscillates; if n<0 such oscillations do not exist 



(see equation (|12|)). To judge whether the oscillating solutions of MDG are physically acceptable, 
one has to find their physical consequences and possible manifestations in the real world. 

2) The dilatonic field <& can change its sign, i.e. in general phase transitions of the Universe 
from gravity (<& > 0) to anty-gravity ($ < 0) and vice-versa are possible in MDG. In general, the 
hyper-surfaces $(a;) = in space-time are dangerous singular Cauchy surfaces: as seen from first of 
equations initial conditions on them do not determine second time-derivatives of metrics and 
evolution. This problem needs careful investigation, but in the case of the RW Universe considered 
here it has a simple positive solution: the second of equations (^) and equation (|^) instead of the 
first of equations (j^) determine the evolution even for initial condition ^{tq) — 0. 

Depending on the form of cosmological factor n(<f>) in some of MDG models the above transition 
from gravity to anty-gravity are possible, but such new phenomena do not exist in models with 
proper cosmological factor. To decide which model describes better the real Universe we first have 
to study theoretically the properties of this rich class of models and then to look for observational 
evidences in favor of some of them. 

3) In the spirit of Max principle Newton's constant depends on the presence of matter: G ~ 
l/^n^l/e, see equation (|l5). 

4) For simple functions a(r) the cosmological factor n($) and the dilatonic potential V{^) may 
show unexpected catastrophic behavior including terms ~ (A^)"^/^ (A$ — $ — $(A*)) in vicinity 
of the critical points A*: ^'(A*) = of the projection of analytical curve {n(A), ^(A), A} on the 
plain {n, The conditions n'(A*) = 0, n"(A*') = 0, n"'(A*) ^ 0; $'(A*) = 0, $"(A*) ^ ensure a 
behavior 

n($) = e"'^^* (n* + 03/2(A$) + 02(A$)) , 

where fT(A) := (3(1 + ^)h^ - e/$) / (joo - 

The potential V{^) has the same critical points $* as the factor n($) and in numerical study of 
above examples shows the following behavior: a local minima around the point 1> and critical points 
at one, or at both ends of some interval [$/e^f , ^rightl ^ ^- interesting open problem is to find 
the general conditions on the scale factor a(T) which make impossible these catastrophes, yielding 
an everywhere analytical cosmological factor n($) and potential V{^) with desired properties. 
The existence of such class of scale factors a(r) is clear: in the direct cosmological problem we 
can choose a desired functions n(<i>) and V{^) and then we will obtain the corresponding function 
air). 

5) Clearly one can construct MDG model of Universe without the typical for GR initial singu- 
larities: a(To) = and with any desired kind of inflation. 

6) Because dilaton field is quite massive, in it will be stored significant amount of energy. In 
RW model this energy has a homogeneous and isotropic distribution in the Universe. An interesting 
open question is: may the field $ play the role of a dark matter, or of a dark energy in the Universe? 

7) A very important problem is to reconstruct the cosmological factor n($) of the real Universe 
using a proper experimental data and the astrophysical observations (see |12| where such problem 
was considered first in the framework of more complicated models with two unknown functions). 

8) As a further steeps for a justification of the parameters of the MDG model and its relevance 
to the physical reality one should study the possible consequences of the existence of the dilaton $ 
for binary pulsars (which give precise test of the extended theories of gravity) and for the curves 
of the rotation of galaxies which are not explained by Newtonian theory of gravity and GR. 

It's clear that MDG offers new curious possibilities which deserve further careful investigation. 
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V. APPENDIX 

Here we give a more extended explanation of the derivation of the basic dynamical equations (|^) 
of the MDG for RW Universe. The first of these equations is a direct representation of the equation 
written for the case of RW metric. The second one represents the (0,0) component of the first 
equation of the system (H) for the same case. The only nontrivial additional equation which one can 
derive from these generalized Einstein equations, considering their (i,i) components, is equivalent 
to the second equation of the system (g). (Note that in the case of RW metric one obtains the 
same equations for i — 1,2,3.) For the RW Universe this equation for the field $ = $(r) acquires 
the form: 

* + 3!0+^ = |(e-3pe). (18) 

From the conservation law Vq Tp = 0, applied for the RW metric, one obtains the same relation as 
in GR: ^ (o-^^) — —3pea^- It makes possible to exclude the pressure from equation ([l^). Now, 



just as in GR, it becomes possible to prove that the equation (18) follows from the basic system of 
dynamical equations (|^): differentiating the second one with respect to the time r, using the first 
one and once more the second one, after some algebra we derive the equation (|l8|). This means 
that any solution of the system (||) satisfies the equation (^8|) and there is no need to take it into 
account, solving dynamical problems for the RW Universe. 

One has to stress that the second order dynamical equations for the RW Universe in MDG: the 
first of the equations (^) and the equation ( |l8|) , are equivalent to the Euler-Lagrange equations for 
the (dimensionless) action: 

a„„, ^ J dr (^-$aa2 - a^^a + fca$ - ia3($n($) + e(a))) . (19) 

It is obtained from the original action (|l|) (per unit volume) substituting in it the RW metric and 
dropping out the corresponding common factor with dimension of action (i.e. the cosmological unit 
for action Ac, introduced in Sec.l). The second of the equations (||), being a first order differential 
equation, represents the corresponding energy integral of this Euler-Lagrange system: 

^2 2r 



etotai = apa + $p<i. - (^-$aa2 - a^^d + fca$ - ia^($n($) + e{a))j = 

+ * {Ml^y + (*n($) + e(a))) EE 0. (20) 



Just as in GR, it must be zero, because of the invariance of MDG under general coordinate 
transformations. In the equation ( |20| ) pa = — 2$ad — and p$ = —a^a are the corresponding 
(dimensionless) canonical momenta for the action ([l9|). 
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